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This is a study of the behavior of wave equations in conformally compactified spacetimes suited 
to the use of computational boundaries beyond Scri+. There light cones may be adjusted for com¬ 
putational convenience and/or Scri-|- may be approximated by a “proto-Scri” spacelike hypersurface 
just outside a de Sitter horizon. One expects a numerical implementation to excise the physically 
unnecessary universe somewhat beyond the outer horizon. As an entry level example I study forms 
of the Maxwell equations and causal relations for an outer boundary in that example. 


I. INTRODUCTION 

The focus of this essay is the matching of a computa¬ 
tional model of a source of gravitational waves (such as a 
compact object binary) to the cosmological environment 
through which any emitted signal travels to a detector. 
The transmission phase is well understood and is essen¬ 
tially the propagation of a linearlized gravitational wave 
through a background cosmology which includes such ef¬ 
fects as a cosmological redshift and gravitational leasing 
due to intervening masses, plus the kinetic effects of the 
detector’s motions. For these purposes the geometrical 
optics approximation is often adequate. Although im¬ 
plementing this transmission model in observatory data 
analysis demands skill, it is simpler than the electromag¬ 
netic case where interactions with dust, gas and plasmas 
en route must also be considered. In this large scale back¬ 
ground the source can be taken to be a point source of lin¬ 
earized gravitational waves having some specific antenna 
patterns and wave forms. The more unsettled question 
then is how to extract from a detailed model of a wave 
generation system the appropriate linearized wave. The 
most common assumption is that the wave generator can 
be modeled as existing in an asymptotically flat space- 
time where waves propagating out to future null infin¬ 
ity (Scri-b or will be described there by linearlized 
gravitation and inserted as the point source is the cos¬ 
mological context. This is clearly justihed by the huge 
difference in scales between the size of the wave generator 
and the radius of curvature of its surrounding spacetime 
(the cosmos or a galaxy or stellar cluster). The ques¬ 
tion then is how to efficiently calculate and extract the 
waveforms at 

In the usual Cauchy evolutions on asymptotically flat 
spacelike hypersurfaces, the waves must be recognized at 
large distances from the source, possibly a long time af¬ 
ter the source motions have ceased violent activity. One 
way to overcome this time lag is through the use of a 
retarded time coordinate. This works well in spheri¬ 
cal symmetry u u, but is poorly defined otherwise. 
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The Pittsburgh group has used such characteristic evo¬ 
lution at larger distances while matching to more com¬ 
mon Cauchy evolution in the interior Q. The use of 
hyperboloidal time slices (some generalization of hyper¬ 
boloids in Minkowski spacetime) appears attractive as a 
smooth transition from conventional time slices in central 
regions to asympotically null spacelike hypersurfaces to¬ 
ward inhnity. Were such a slicing to be combined with a 
straightforward Cauchy evolution, the Courant time step 
condition would defeat the advantages of approximating 
retarded time, as the nearly infinite speed of light toward 
infinity would require inhnitesimal time steps. Thus hy¬ 
perboloidal slicings effectively require a conformal com- 
pactification which makes the coordinate speed of light 
finite again Q. 

While many other approaches iii reformulate the 
Einstein equations, I conjecture that the ordinary Ein¬ 
stein equations might be used for a conformally compact- 
ihed metric g^^i, = where g^y is the compactified 

metric and g^v is the physical metric. If the compactifica- 
tion is successful both metrics will have nonsingular Ein¬ 
stein tensors where 17 = 0; the physical metric because 
it has a vanishing Einstein tensor, and the compactified 
metric because the null hypersurface where S7 = 0 is by 
assumption smoothly imbedded in the conformally regu¬ 
lated spacetime with metric g^i,. Thus the Einstein equa¬ 
tions should reduce to the Einstein tensor of g^y set equal 
to an effective stress energy tensor computed from the 
conformal factor S7 as described in 0, Equation 11.1.16]. 
Although these terms contain fl In denominators, this 
effective stress-energy tensor must be regular since the 
conformally regulated metric and Its Einstein tensor are. 
Before studying this conjecture in constructive detail, I 
want to conhrm that wave equations, evolved numeri¬ 
cally by conventional Cauchy evolution schemes, do not 
fail for some unanticipated cause. Thus this paper formu¬ 
lates a problem for numerical study in which the Maxwell 
equations are to be evolved in a conformally compactified 
spacetime that is nearly Minkowskian. 

In addition to the use of hyperboloidal time slices 
and conformal compactlhcatlon, I suggest two additional 
tools to ameliorate the extraction of outgoing waves in 
numerical evolutions. One is the use of artificial cosmol¬ 
ogy to convert the null hypersurface into a spacelike 
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hypersurface in the future of the initial Cauchy hypersur¬ 
face as proposed in Q . This procedure fattens into a 
shell of spacelike hypersurfaces by changing the value of 
the cosmological constant plausible for the present uni¬ 
verse by many orders of magnitude so that the de Sit¬ 
ter horizon is located, not at cosmological distances, but 
merely well beyond the region where the wave generat¬ 
ing mechanism can be influenced by it. In this sense it 
is similar to the computational tool of artificial viscos¬ 
ity which fattens a shock from its physical thickness on 
the order of the molecular mean free path into a much 
larger length which is still small in the range of scales 
that otherwise are met in the physical problem at hand. 
The second additional tool is the tilting of the null cones 
governing wave propagation in ways that ease the compu¬ 
tational behavior in the unphysical (‘Oz’) region beyond 
the de Sitter horizon and beyond . 
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II. BACKGROUND SPACETIME 

As in the background metric for this study will be 
the de Sitter spacetime 

= -dT^+dX'^+dY^+dZ'^ + {R^/L^){dT-dRf ( 1 ) 
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with an artificial cosmological constant A = -|-3/L^. We 
are interested in the limit 1/Lp = 0 when this becomes 
Minkowski spacetime, and the cosmological constant is 
used to make boundary conditions numerically simpler 
at the de Sitter horizon than they would be at flat space¬ 
time’s . The de Sitter horizon is the null hypersurface 
i? = L in this metric. 

As before I introduce the coordinate changes 


T 

— =u + 
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(2a) 


yield the regulated de Sitter metric, but could be mod¬ 
ified outside r = 2 to control error propagation in the 
unphysical region. 

To see the behavior of the light cones for this regulated 
metric, we calculate the coordinate speed of light dr/du 
in radial null directions. This gives 

_r /L'^)S ±ao 

1 + (s2/l2)52 • 

For the de Sitter values above this reduces to 


s 
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and 


■Cout = (1 + ^r)^ 


( 8 ) 


where R^ = X'^+Y'^+Z'^ = and = x'^+y'^+z'^ = 
x'a;*. This brings in to r = 2 in the Minkowski case 
and in the de Sitter case makes this a spacelike hypersur¬ 
face beyond the de Sitter horizon which I call proto-.y'*'. 
The hypersurfaces of constant u are then hyperboloids 

[T - s{u - l)]^ - R^ = (3) 

in Minkowski spacetime and are also asymptotically null 
spacelike hypersurfaces in the de Sitter modification. 
This leads to a metric which is singular at r = 2, but 
only in a conformal factor s^/q^ = s^/(l — which 

does not appear in the Maxwell equations (see |3|). Thus 
we can drop this conformal factor and our test problem 
is to solve the Maxwell equations in the resultant metric 
which is of the usual 3-1-1 form 

ds^ = —Q^dt^ + 7y (dx* -I- P^dt){dx^ + P^dt) (4) 


-(1- l4.)2 + (,2/j^2)(i+l^)2^2 

1 -f (s2/L2)5'2 


( 9 ) 


Note that the constant r hypersurfaces are never space¬ 
like when s/L = 0 since then the inward speed of light 
is Un < 0. When the de Sitter option is used one finds 
Vin > 0, i.e., both the inward and the outward sides of 
the lightcone point toward increasing r, for an interval 
around r = 2. 


III. TILT 

In evolution codes which permit setting the computa¬ 
tional boundar y o n a coordinate sphere (such as those 
described in [ifl[illUS ) ^^^e outer boundary could be 
set at r = 2 or slightly beyond in the metric described 
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above (or ones asymptotically similar). However, if a cu¬ 
bic boundary enclosing is to be used, there are addi¬ 
tional behaviors beyond that need attention. For ex¬ 
ample, if the artificial de Sitter horizon is set at L = 10s, 
then the r = constant hypersurfaces are spacelike only 
in the narrow interval of (approximately) r = 1.81 to 
r = 2.21. But a boundary cube enclosing r = 2 has cor¬ 
ners at r = 2-\/3 « 3.46. Thus on most of such a cubic 
boundary the light cones (from equation would allow 
inwardly propagating waves, requiring careful attention 
to the boundary conditions which one hopes to avoid. 
In this case I suggest that the metric be modified out¬ 
side r = 2 so that all r = constant hypersurfaces beyond 
r = 2 are spacelike — i.e., so that Um > 0 in the en¬ 
tire computational domain beyond r = 2. Analytically 
this would have no effect on the physical solution inside 
r = 2 provided all fields propagate causally. It is also 
straightforward to do in the test case of a Maxwell field 
with a fixed background metric. If these methods were to 
be developed for the Einstein equations (without using 
a spherical boundary such as r = 2) it would necessitate 
distinguishing, in the Einstein equations beyond r = 2, 
the metric begin evolved from the appearances of the 
metric as coefficients of the principle derivatives which 
determine the causal relations in the evolution algorithm. 

In the Minkowski case s/L = 0 (no artificial cosmol¬ 
ogy) it is not possible to make a change in the metric 
only beyond r = 2 which avoids inward pointing null rays 
just beyond r = 2 — in that case one has, from the r < 2 
formula above for Vin that {(P/dr‘^)vin = —^ at r = 2 so 
that fin must have a local maximum value of zero at that 
point and is necessarily negative (inward) in a neighbor¬ 
hood of that point, including some r > 2 points. Eor 
the de Sitter case there is again a maximum of Uin near 
r = 2, but its value there is positive so it can be smoothly 
turned upward (to avoid zero or negative values) in suit¬ 
able intervals beyond the de Sitter horizon (and beyond 
r = 2 for small values of s/L). 

IV. MAXWELL’S EQUATIONS WITH 
CONSERVED CONSTRAINTS 

We take the Maxwell equations in the form given in 
the DeserFest paper [^. The quantities involved are all 
various components of the usual 4-dimensional Maxwell 
fields F^i, and In particular we 

define, S'”* = = [ijk]Hk and = Ft, 

Fij = Bij = . Here [ijk] is the completely anti¬ 

symmetric symbol \ijk] = 0, ±I with [123] = -l-I. The 
index 0 here refers to the coordinate basis vectors dt or 
dtt. The fundamental fields in the formulation are IB* and 
J)*. The constraints read 

= 0 = . ( 10 ) 

In terms of these primary fields, two auxiliary fields Ei 
and Hi are computed by the formulae 

E^ = (a/V^) 7 yD^' + [ijk]l3m'^ (Ila) 


and 


(11b) 

The evolution equations are then 



9t$® = -[i3k]d,Ek 

(I2a) 

and 

dt'S^ = [ijk]d,Hu . 

(12b) 


Because the combinations of metric components appear¬ 
ing in equations OH are each conformally invariant, it 
follows that field components IB®, J)®, and Hi satis¬ 
fying these Maxwell equations are also each conformally 
invariant. 

By differentiating equations 03 one find equations for 
the evolution of the constraints 

9t(5,D®) = 0 = 5t(a,»®) . (13) 

These show that the constraints are conserved, so that if 
they are initially zero they remain zero within the limits 
of the numerical approximations used. This sort of be¬ 
havior has been found less than optimal as shown, e.g., 
in a Maxwell example by Baumgarte et al. An alter¬ 
native form for the equations is given in section^T] below 
which emphasizes use of the direction normal to the time 
slices. 

V. ENERGY CONSERVATION 

Because the metrics here are stationary with Killing 
vectors ^ = do = dt there is a conserved energy. The co¬ 
variant “conservation” law = 0 leads to a true con¬ 
servation law dvp'' = 0 where p®® = which 

gives a conserved integral / p* (Px whose changes can be 
monitored by the energy flow through a boundary inte¬ 
gral. Most of these equations hold only in a coordinate 
basis, so we express the energy current components us¬ 
ing the coordinate basis components of the Maxwell field 
from section nvi above. 

For the electromagnetic field the stress-energy tensor 
density is given by 

- IS/F^pr^ ■ (14) 

When the components in 3-1-1 form are inserted with = 


(I; 0,0, 0) one finds for the energy density 


4V = i(S)'=Efe + <B'®i7fc) 

(15) 

and for the energy flux 


47rp® = [ijk]EjHk . 

(16) 


These expressions, in addition to the constraints, pro¬ 
vide possible checks to monitor the validity of numerical 
implementations of the Maxwell field evolution. Note 
that these expressions are conformally invariant (since 
each field is) and yield spatial coordinate independent 
integrals as no further metric factors enter when forming 
integrals of these densities. 
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VI. MAXWELL’S EQUATIONS WITH 
DISPERSING CONSTRAINTS 

I here reformulate the Maxwell equations from the 
DeserFest paper and section nYi above to make the 
principal terms invoke the derivative along the normal 
Qn = dt — P^dk- This makes it natural to work with 
components of the fields taken in a frame whose time 
axis is normal to the t = const hypersurfaces. Thus I use 
the basis 


and similarly 55^ = ^[kij]Bij. 

As the fundamental fields being evolved remain J)® and 
Bij or 18^, the constraints are the unchanged equations 
m above. Evolution equations in terms of the frame 
based components are found by substituting equations 
(im into the evolution equations OS- 

In the slicing normal frame 113 it is straightforward to 
relate the ^ components to those of F since there are no 
time-space metric components gNi in that basis. Thus 
we find 


CN = On = dt- 13'" dk , (17a) 

e^ = d^ , 

and dual basis 

= dt , (17b) 

oj>^ = dx'"+j3'"dt . 

Note that this basis is constructed in a conformally in¬ 
variant way since the shift /3® is a conformally invariant 
aspect of the metric tensor. Tensor components will be 
identified by a tilde when referred to this basis. Thus 
the metric components g^u can be read from equation 
0 and are gNN = -o?, gN% = 0, and gtj = = gtj. 

For the other fields of interest we define 

A = —tjjdt -I- Aidx" = -I- AiLo" , (18a) 


-Eiu;^ Alu'f A uj^ , 

(18b) 

5" = S'cat a Cj -I- \h^^ei A Bj . 

(18c) 


By inserting the basis vectors from equations m these 
can be expressed in terms of the coordinate based com¬ 
ponents used in @ and section Hvl and one finds that the 


principle fields are the same: At = At, Bij = Bij, and 
T)® = Db The subsidiary field components do change 

= f3'"Ak , (19a) 

Ei = E,+B,jP^ =E, + [ijk]P^^'" , (19b) 

and 

or 

H, = H,-[ijk]f3^^'" . (19d) 


The antisymmetric 3-tensors have been renumbered in a 
metric independent way to 3-vectors of a different density 
by relations such as 

Hk = ^[kiM^ or S)"^ = [ijk]Hk (20) 


Bzj = (l/av^)7»fc7ii'ft''* (21) 

or 

Hi = {a/y/^)jij^^ . (22a) 

and 

Ei = {a/./^)^ij& (22b) 

To obtain equation 122all from ea does, however, need 
the definition of a determinant in the form 

[i3k]l^ploqlkr = [pqr\l ■ (23) 

The evolution equations for ^ and S) found in this way 
are 

{dt - £/ 3 )»* = -[ijk]djEk - p^dk^’") (24a) 

and 

{dt-£0)^" = [ijk]d,Hk-P"{dk^'") . (24b) 

In these equations the Lie derivative can be expanded as 


= dj{P^^P-^^djP' (25) 

and similarly for These equation are un¬ 

changed, just rearranged, from OS- But on the right 
hand side of each of equations the terms containing 
P vanish in view of the constraints. Thus these terms can 
be omitted, leading to 

{dt - =-[ijk]djEk (26a) 

and 

{dt - £^)S)* = [ijk]d,Hk . (26b) 

This system (I2SI is physically equivalent to 113 - i.e., 
agrees with them when the constraints are satisfied — 
but may have different behavior in numerical approxi¬ 
mations when the constraints will not be exactly zero. 
This is in fact the case, as the evolution equations m 
leave the constraint values constant in time, while the 
evolutions system show a different behavior as I now 
show. 
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By taking the divergence of equations IIH) one finds 


{dt-p^apidk^^) 


(27a) 

(5t-/3*5,)(5fcD'=) 


(27b) 


These are each ordinary differential equations for the re¬ 
spective constraints along curves with tangents {dt—P^di) 
along which the constraints will have an exponential be¬ 
havior with e-folding time (9^/3*)“^. Although all the 
equations in sections CYI El and ED are valid for any 
metric of the form 0, we are most interested in the reg¬ 
ulated de Sitter metric defined by equations 0 through 
There one calculates that (for small values of s/L) 
to good approximation /3® = —x* and dkP^ = —3 so that 


the constraint propagation equations m show the con¬ 
straints propagating along the normals to the constant 
u hypersurfaces while decaying exponentially with an e- 
folding time of Au = 1/3. It may not be possible to 
preserve this favorable constraint behavior if the metric 
is modified by tilting the light cones beyond . 
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